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Abstract 

A recently introduced Renormalization Group approach to frustrated spin models is applied in 
three dimensions through Monte Carlo computations. A class of spin glass models is analysed, 
with correlated disorder variables given by a Z2 gauge field. Evidence is provided for the 
influence of deconfinement phase transition of gauge fields on the behaviour of the associated 
spin models, namely in the transition from glassy to ferromagnetic behaviour. Universality 
classes are determined on fundamental ground by characterizing the fixed point. Moreover, the 
RG analysis provides good estimates of the critical temperature and the thermal index v, with 
moderate computer time. 
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group methods. 



1 Introduction 

The renormalization group (RG) is an invaluable tool in the understanding of critical phenomena [|IJ |2||. 
It provides accurate estimates of critical indexes as well as a theoretical foundation of the universality 
principle and of the scaling theory. Many authors |l3|| have profitably used these ideas in spin glasses 

(SG's), even if a general RG framework for SG's was incomplete. 

In a previous paper [|l8]] a RG approach suitable for SG's has been proposed and fruitfully exploited 
through Monte Carlo (MC) computations in dimensions d = 4. A sketch of expected universality classes 
and the role of frustration were discussed. In the present work, these ideas has been fully developed in the 
more relevant and harder case of d = 3. The coarse graining transformation is carried out on the overlap 
field. Because of disorder average, overlap probability measure does not have the usual form exp(— H) , 
hence it is not a priori clear which kind of additional interactions are generated by the RG transformation. 
Our proposal is to keep Hamiltonian fixed while adding gauge invariant terms to the disorder distribution, 
so that spins get effective correlations via correlations of gauge fields. Asymptotic decrease regimes of 
Wilson loops in pure gauge model [|J, namely the weak and strong decrease, can induce very different 
behaviour on related spin systems. In the former case, frustration effects are very small locally and increase 
weakly at long range: spins prefer to order ferromagnetically. In the case of strong decrease, instead, 
average Wilson loop rapidly fall to zero and a SG behaviour is expected. The existence of two fixed points 
corresponding to the above mentioned universality classes is supported by evidence from our analysis. 
Moreover, our estimates of critical temperature T c and thermal index v are in good agreement with previous 
results [ |i"4] , |l7[ ]. Notice that a less frustrated system gets ordered more easily, giving rise to higher critical 
temperatures. This property can be exploited to perform Monte Carlo computations of SG models at very 
low T/T c ratios |p|. 
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Numerical methods based on Monte Carlo simulations of SG's are particularly hard in d = 3 basically 
because of the low temperatures involved. This requires very long Monte Carlo runs to reach ergodic 
regime. Parallel Tempering (PT) algorithm, which we implement, is adequate to face this kind of problems 
but unfortunately it hides the meaning of MC autocorrelation time, making standard error estimates not 
efficient. Aiming to control both ergodicity and MC-errors, we employ a method based on the analysis of 
MC-PT dynamics, which turns out to be quantitative and more stringent than usual PT efficiency checks. 



2 The spin glass model and the overlap field 

Let <7 X be Ising spins located at the sites of a d-dimensional cubic lattice A, with L points on each side 
(x € A C Z d with toroidal topology, i.e., periodic boundary conditions are assumed). The typical spin 
glass model of Edwards and Anderson [Q] is defined by two prescriptions: 1. The spins are distributed 
according to the Boltzmann weight given by the following Hamiltonian 

H(J, a) = - J xy a x a y ; (1) 

2. The quenched disordered interactions J xy are random variables, with a given distribution. The sum in 
(|l|) is over pairs of neighboring sites in A. We will denote by E the expectation on J variables, and define 

E(J xy ) = 0, E(J 2 xy ) = l. (2) 

The Boltzmann-Gibbs measure on the spin variables will be denoted by angular brackets (•). Averages 
over the disorder are taken only after Boltzmann averages are calculated, and the thermodynamic limit for 
the appropriate quantities is eventually taken afterwords. A generic observable can be expressed in terms 
of correlations of correlation functions, 

E[(h(a))... (f .(*))], (3) 

where fi, i — 1, . . . , s are functions of the spin field. Consider the following gauge symmetry of the 
Hamiltonian, 

Jxy — * J X y = £x-Jxy£y ^ 
a x — >• a' x = £xO~x 

where e x — ±1 are the gauge group parameters. If the distribution of the J's is gauge-invariant, the model, 
characterized by observables (^), is also gauge-symmetric. This happens, for instance, in the traditional 
case of even, independent, identically distributed J variables (E-A model). 

A central concept is the overlap field. Consider s replicas (copies) of the spin variables 
1, . . . , s. The Hamiltonian of the replicated system is given by H = ^2 a H(J, a^): replicas are indepen- 
dent from each other, but feel the same disorder configuration J. For any pair of replicas (a,b), 1 < a < 
b < s, define the overlap field, q x a,b \ 

q x a ^ = e Z 2 . (5) 

The probability distribution of the overlap fields, /i, can be implicitly defined through the overlap expecta- 
tions, that involve both the thermal average and the average E over disorder. For any smooth function F 
define 

(F (q^\q y ^\. ..)))= E [(F 4 2 >4 3 \ ...))], (6) 

where expectation with respect to the /i distribution is denoted by ((■)). All physical observables can be 
expressed in terms of overlap observables, so that the full physical meaning of these models is contained in 
the overlap probability measure. 



2 



Let us introduce the volume average of the overlap field, the total overlap, often simply referred to as 
"overlap", 

9 (Q ' b) = ^£^' b) - ^ 

For two replicas, the distribution of the overlap will be denoted by . An other interesting observ- 
able is the two points, connected correlation function of the overlap field: 

' ' X 

3 A renormalization group for spin glass models 

The total overlap is the order parameter of the model. In the high temperature phase, including the critical 
point, it should be zero in the infinite volume limit, while in the low temperature phase it is expected to 
fluctuate |Q, ^, [uj. Therefore, the critical point can be characterized by a divergence of the correlation 
length associated to (|^). For this reason, it is natural to define the RG transformation on the overlap field. 
In the following we consider only two replicas, and omit the replica indexes on the overlap variables. 

Let B™ C A be the cube of size n located in x and A n C Z d a cubic lattice of side L /n. Because of 
Z2 symmetry, the coarse grained field q' x , can be defined [||] on A„ as the sign of the majority of overlaps 
inside a block centered on x = nx'. The renormalization group transformation for a rescaling factor n then 
acts on the overlap distribution /i as follows: 

Atf) =5>(9) 11 8 I q' x , ,sign I }_J ,/,,, j j . (9) 

{<?} x'eA n \ \y£ B Z 

where the sum on the r.h.s. runs over all the 2' A ' configurations of overlaps q x and S(a, b) is 1 for a = b 
and zero otherwise. For m — n k , one defines the distribution {q^) as the result of k iterations of (^]), 
and the semi-group property holds by construction. 

Extending these definitions, in the appropriate way, to the thermodynamic limit, the transformation may 
be iterated indefinitely and eventually a non-trivial asymptotic distribution /1* may be reached. However, to 
carry out useful RG calculations an approximation scheme must be chosen. In this work we have used the 
Monte Carlo approach, for which the main approximation is the use of a finite lattice. The basic effect is that 
of neglecting the influence of distant regions of the lattice on each other in the calculation of /J,' n {q') from 
n(q) [||]. Indeed, a fundamental property of a RG transformation is that the parameters representing the 
renormalized distribution must depend analytically on the unrenormalized parameters. One assumes that 
the dependence of p! n on /x is analytic, because local observables of q' in one region of the lattice are not 
appreciably affected by distant regions of the system, in spite of the eventually infinite correlation length. 
For the very same reason, the computation of local observables of the field q' should be weakly affected by 
finite lattice effects, on the contrary of what happens to observables of the unrenormalized field q. 

On a finite lattice the transformation ([))) may be limited to a very small number of iterations. Another 
typical problem is the truncation error due to the projection of renormalized distributions onto a space with a 
restricted number of couplings. To minimize the effects of these approximations, a different RG scheme [|fU 
|l8| ] has been adopted. We define the renormalized fields q' x on A„ as given just by the transformation (g) 
with n = L/l (with small typical values, e.g. t = 2, 3); the renormalized distribution will be denoted by 
n'(q'). Notice that, with only one iteration, A„ is a lattice with volume £ 3 . Due to the smallness of A„ only 
a few couplings are allowed and there is no truncation error. The procedure is repeated for different values 
L, L' and the ratio L/L' is interpreted as a scaling factor [|llj (but see also [|i"2|]). Moreover, the sequence 
will give better results for increasing n, because the group parameter will be large and, at the same time, 
the size of the initial system will be affected by smaller finite-size effects. Given the large size of blocks, 
we also expect the results to be not dramatically dependent on the choice of blocking transformation. 
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The RG transformation is naturally defined on the probability distribution of the overlap field. For the 
purpose of calculations this distribution is well characterized by the expectations of a set of (translation- 
invariant) observables. It is the approach that has been followed in this work. However, the question may 
be raised of how to express such a distribution in terms of physically meaningful parameters. It is clear that 
the traditional exponential form, exp(— H), has no physical relevance, due to the involved definition of the 
overlap field distribution in the original model. To state the question in a different way: which additional 
interactions between the microscopic variables are generated as a result of the RG transformation ? The 
answer to this question is important in order to get high precision estimates of the critical indexes, through 
the observation of the effects of irrelevant perturbations to the fixed point. We propose the following pa- 
rameterization, that proved to be effective in the case of dimension d = 4 [jig]. The Hamiltonian, given by 
([]]), is kept fixed, while the disorder is distributed according to a general gauge action: 

PK (J) = exp (j2i Ki W ^) ' (10) 

where the K t EM. are parameters and the Wis are Wilson's loops, i.e., products of J's along closed paths. 
Let us stress that the effects of this parameterization are that spins get non trivial extra correlations through 
the "dressing" of interactions. The resulting set of spin models will be discussed in more detail in next 
section of this paper. 



4 Spin models with quenched gauge-field interactions 

Consider the E-A spin Hamiltonian ([!]), and assume the quenched interactions are distributed according to 
a general, gauge-invariant distribution function such as (|1C|). A fairly large set of models may be defined in 
this way. For the sake of simplicity, let us restrict our attention to the disorder distributions of the following 
form, 

Pk{J) = C K exp K x J ly + K * E 4» + K * 12 D <* I < ^ 

\ (x,v) (x,y) a J 

where Ck € M is a normalization constant, and the symbol □ denotes the plaquette terms of the kind 
Jx,yJy,zJz.wJw,x- Expectation with respect to this distribution will be denoted by Ek- We notice that 
the Gaussian Edwards-Anderson model corresponds to K\ = —1/2, = K3 = 0, while the Z2, E-A 
model is obtained in the limit K\, K2 00, s.t. K1/K2 = —2, K 3 = 0. But a more interesting model is 
obtained by adding a plaquette term (i.e., K 3 ^ 0): 

p K {J) = C K e K ^« a ^ J€Z d 2 lAl , (12) 

Considering the disorder variables only, this is the well known pure gauge Z2 model [S], The partition 
function, the free energy and the internal energy density are defined by 

Z A (K 3 ) = E { JxH=±1} exp (K 3 £ Q a a ) = e -K,F A (K,) (13) 
u A (K 3 ) = ~E K3 (E q n a /N n ) ■ (14) 

The energy density u\ is equal to minus the average plaquette, thus it gives a measure of the short - range 
frustration of the spin system. 

Let us denote by w 1 the product of the J's along a closed path 7, 

w 7=n j ^"ii n «' (i5) 

7 aes 

where S is a surface bounded by 7. Two canonical asymptotic regimes may be distinguished: a weak 
decrease regime, characterized by E(u>j) ~ e~ L ; a strong decrease regime, where E(u>j) ~ e~ A . L and 
A denote the perimeter of 7 and the area of the surface S. In various models a transition from the strong to 
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the weak decrease is observed as the average plaquette increases. The model (|l2|), in 3 dimensions, exhibits 
a second order phase transition at A3 — A3 ~ 0.7613. The corresponding order parameter is the square 
Polyakov loop, p 2 = (p\ + p\ + p 2 )/3, where pi is the average Wilson loop along a path that closes itself 
exploiting the periodicity of the b.c.'s in direction i. The quantity E(w 7 ) is relevant for the spin system 
associated to the gauge field [^], as it is related to the average frustration on the range 7. 

We are now in the position to sketch the behavior of the spin system corresponding to the pure-gauge 
disorder distribution ([jj]) [Q: 

A3 — > 00 The disorder is in a ground state of the gauge Hamiltonian, characterized by Ek 3 (0) = 1, 
Pi = ±1. Up to a gauge transformation, the spin model is an Ising model, with periodic b.c. in the 
directions i for which p t = 1 and anti-periodic in the others. 

A c < A3 < +00 The disorder is in the weak decrease phase: frustration effects are very small locally 
(Ek 3 (0) ss 1) and increase weakly at long range. The spin system is expected to be ferromagnetic. 

— K c < A3 < K c The disorder is in the strong decrease phase. The average plaquette is approximated by 
Ek 3 {0) « A3, except in the region very close to the transition. However, frustration moves rapidly 
toward 1/2 for increasing range. Spin glass behavior is expected. 

A3 — > — 00 The gauge field is in a "fully frustrated" configuration, i.e., all plaquettes are frustrated. 

The behavior in the region — 00 < A3 < K c is still an open question. 

These different qualitative regimes should correspond to different universality classes of the spin sys- 
tem. Exploiting the RG approach introduced in the previous section, we will be able to characterize the 
universality classes in terms of different fixed points. 

5 Analysis of Monte Carlo computations 

The Renormalization Group analysis described in section [}] was performed numerically by dynamic Monte 
Carlo computations. We employed Hamiltonian ( pj) in dimension d = 3, with periodic boundary conditions 
and Z2 quenched disorder distributed according (|12|). Besides the traditional Edwards-Anderson model, 
corresponding to the case A3 — 0, we simulated two other models, with A3 = 0.3 and 0.8, in order 
to verify the conjecture on universality classes made in section ^. Moreover, the MC-RG analysis was 
performed on the Ising model, corresponding to A3 = 00, as a test of the method, and for a comparison to 
the other models. 

The choice of A'3 values deserves some comments. A3 = 0.3 is far from the critical value Af ~ 0.7613 
to avoid cross-over effects, and is quite different from the E-A model: frustration is /0.3 ~ 0.35, compared 
to fo = 0.5. In the deconfined phase, A3 = 0.8 gives a very low frustration /o.g — 0.015. Lower values of 
A3 provide slightly higher frustration, but we considered unwise a value too close to A3, even though the 
critical coupling is lower in a finite lattice. Nevertheless, the model is quite different from the Ising model, 
as is shown, for instance, by the value of the square Polyakov loop, p 2 ~ 0.65. Among the eight possible 
pure phases, we chose the one characterized by positive Polyakov loops, i.e., pi > 0. 

For each model, a quite large number of samples (with different configurations of the disorder J) were 
simulated, as reported in table |l[ The required disorder configurations were generated by an independent 
MC process on the gauge fields, according to distribution (|T2"|). This step is quite simple with nowadays 
computers. The average correlation between J configurations is never greater than 10 -4 . 

The MC runs on the spin variables are more complicated. We simulated lattices of linear size L = 
8, 12, 16 for all models, together with additional L = 24 for models A3 = 0.8 and Ising and L = 
10, 14 for A3 = whose significance will be clear later. The code makes use of Parallel Tempering 
(PT) algorithm and multi-spin encoding. For every Metropolis update, a PT step was performed. Two 
replicas were simulated through independent MC chains in order to compute overlap observables. The 
RG transformation was carried out following the procedure explained in section || JTT|]. Therefore, the 
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RG parameter n takes the values 4, 6 and 8 (and 12) as well as 5, 7 for L = 10, 14. To characterize the 
distribution //, we measured the following observables of the random field q': 



1 



E MM 



A 4 =(S 2 )) 



(16) 



3|A'| 



A 2 



1 



E fe^x+i) 



^>=«s 4 » 



(17) 



3|A'| 



x,i:|i|=A/2 



^3 



1 



E (ixlylWy 



» 



A 6 =«(0.4 + 5 2 )- 1 )) 



3|A'| 



(18) 



where S — |A'| _1 Qx> an( l tne sum i n ^3 mns over tne 4-uples of sites located on plaquettes. Notice 
that observable Aq is sensible to the small values of S, as opposed to A4 and A$, and 0.4 is an arbitrary 
cut-off. For the sign function involved in the definition of q', a tie-breaker was used in the undetermined 
zero cases. Models with L = 10, 14 are unaffected by this procedure and, in principle, this difference in the 
transformation rule could provide different results. Moreover, we have explicitly measured the derivative 
of each observable, computing the connected correlation function with the Hamiltonian. In the definition 
of the blocks, the choice of the origin is arbitrary. In order to improve the quality of statistics, we choose to 
average the measures over eight different origins. The lower simulated temperature, T m j n in Table [j], was 
chosen around 0.9T C , according to previous estimates of T c , for the K3 = model. The efficiency of the 
PT algorithm was checked through the acceptance rate of temperature swaps. We ensured it was roughly 
1/2 in the most difficult runs, and never below 40%. 

In the rest of this section we discuss the checks on ergodicity of MC runs and error estimates. Generally 
speaking, as higher values of K3 correspond to lower frustration, systems should get ordered more easily 
so that we expect critical temperatures to increase with K3, carrying on a global improvement of MC-PT 
performances. On the opposite, negative values of K3 presumably give rise to very low critical temperatures 
thus making simulations very difficult, at least in dimension d = 3. A well known problem in performing 
Monte Carlo computations on glassy systems is the difficulty in reaching ergodic regime, due to strong 
metastability effects and to the large number of nearly degenerate ground states. A very powerful — albeit 
indirect — criterion to ensure ergodicity makes use of the observation of (Monte Carlo) time spent by the 
system at each temperature of the PT set. This must be approximately the same for each temperature. We 
measured the ratio r of the time spent in the most visited T over the less visited one. The histograms of the 
results are presented in Table [j]. 

Ergodicity may be verified explicitly by analyzing the time sequence of the measures. By the same 
token we can estimate the statistical error of the measures taken on a single sample. The algorithm employed 
measures observables at each MC-PT step. These measures were averaged in groups of 1024, providing a 
set of a few hundred points per sample and per temperature. One problem of PT algorithm is that of properly 
defining an autocorrelation time. However, these points may be considered reasonably independent from 
each other. Thus, we can estimate the statistical error on a single sample. Then the procedure is repeated: 
the measures are averaged over groups of points, obtaining a different estimate of the error. Clearly, as the 
size of the groups increases, the resulting points will be less correlated producing a better (usually larger) 
estimate of the error. In Figure [l] (left), we plotted the histogram of the single-sample relative errors on 
observable A2, in our worst-case simulation (if 3 = 0, L = 16, T = 1.02). The quantity I denotes the 
number of initial points and a is the number of groups so that the size of groups is roughly I /a. The results 
seem to be very stable in a. 

We can check ergodicity in a similar way. On each group i = 1, . . . , a we compute, along with the 
average observable xf\ the error on the average, i.e., its standard deviation (denoted by t^). If ergodicity 
is reached on a time scale I /a then the x\ s will be nearly aligned and fluctuations around the mean will 
be of order ~ f = a -1 ^ t^. On the other hand, if //a is a too short time scale, the points will spread 
with very large fluctuations around the mean. Let us denote by a a the standard deviation of group averages 
xf\ To characterize ergodicity we thus measure the ratio r/ = <t /t, expecting r) ~ 1 when ergodicity is 
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Table 1: Parameters of MC-PT runs. The same number of MC-PT steps, -/Vmc-pt, was used for thermal- 
ization and measurements. We report also the number of different disorder samples Ng, the total number of 
temperatures Np ( together with minimum and maximum values) for the PT algorithm, and the histograms 
on visiting frequencies. 

reached, and 77 3> 1 when is not. The histograms of the results for different samples are shown in Figure [j] 
(right). Let us underline that the results should be extrapolated to a = 1: a "good" MC run should not be 
much longer than necessary. 

From general arguments, the total statistical error on observable A, AA, can be estimated as follows: 

(AAf = $t±?- (19) 

where e| is the sample-to-sample variance, N§ the total number of simulated samples and a is the single- 
sample error. As typical sample-to-sample fluctuations of observables are of order 100%, the single-sample 
errors can be neglected. 



6 Numerical results 

Performing the Renormalization Group analysis explained in section ||, we are able to characterize uni- 
versality classes depending on frustration parameter, pointing out the influence of deconfinement phase 
transition of gauge fields on the behavior of the associated spin model, namely in the transition from glassy 
to ferromagnetic behavior. In particular the results confirm the conjecture [||| that the asymptotic regime 
of long range frustration is responsible for the glassy behaviour. Universality classes were determined by 
characterizing the fixed point, therefore on fundamental grounds, without relying on phenomenological 
comparison of critical indexes. Moreover, the RG analysis provided good estimates of the critical temper- 
ature, T c , and thermal index v, of the four models considered: K3 = 0, 0.3, 0.8 and Ising. As we will 
show in this section, the fundamental determination of universality classes allows to consider and 1/0.3 
as independent estimates of the thermal index of the glassy fixed point. 

In order to damp irrelevant perturbations, and to enhance the effect of the relevant coupling [pT[, we 
chose to employ the largest possible values of n = L/L' = 4, 6, 8, (12), thus setting the renormalized 
lattice size equal to L' = 2. As a drawback, observables A lt . . . , A 6 turned out to be strongly correlated, 
making it impossible to perform a full linearization of the RG transformation on the fixed point and to 
estimate irrelevant exponents. 

Observables defined on different lattices characterize the trajectory of RG flow and, at the critical 
temperature, they move toward a fixed point where they become independent of the initial lattice. As a 
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Figure 1: Relative error (on the left) and ergodicity distribution (right) of observable Ai on various time 
scales 1/ a in the E-A model K% = 0, L = 16, T = 1.02. The total number of samples examined is 1024. 
Relative error diagram is quite stable and more than 99% of sample have relative errors lower than 0.08, 
on all time scales. On the right, the distribution of rj shows that on short time scales (a — 26), 63% of 
samples have r\ > 2; this number drops to 41% for a = 8, and to 25% for a — A, indicating run length is 
essentially correct. Single sample errors and ergodicity both improve at higher temperatures. 




function of temperature, observables will cross in a point which can be estimated as the critical temperature 
T c . The linearization around the crossing point of two quantities A(L\) and A(L,2) determines the thermal 
exponent through: 

\ dA{L 1 ) I / dA(L 2 ) I 
dp l/Sc / dp 1/ 



In 

-\A) = - 



HL1/L2) (20) 

where /3 _1 — T and derivatives are computed at the critical point. In Figures [| and |3| is plotted observable 
Ai as a function of temperature, the points size being larger than the error bars, except in the detail of 
Figure ^| In accordance with Jl4||, Figure ^ seems to rule out a Kosterlitz-Thouless transition type, for in 
that case the curves should merge under T c . Moreover, the resulting estimates of v obtained with formula 
( |2C| ) are the best proof that curves cross with different slopes. 

Data, for each Ai, were fitted to allow for extrapolation to the critical point. We notice incidentally that 
the measures of derivatives with respect to greatly improved the precision of the results. The errors on T c 
and v were estimated with the following procedure: measures taken on different samples were grouped in 
32 blocks and averaged over 31 blocks to compute the results. The procedure was repeated with different 
blocks, obtaining 32 different estimates of T c and v. Then, errors were estimated with jackknife formula. 

The results shown in Table || were averaged over the six observables. Notice that, for glassy models, 
only two of the three estimates provided are independent, as they all come from the same three lattices. 
The error on T c and, as a consequence, on v is larger for the pair (16, 12) because of the very close slopes 
of curves. A transient effect due to the flow of irrelevant couplings seems to be present, but cannot be 
measured within current precision and lattice sizes. Results for K3 = 0.8 have much higher precision, 
revealing that sample-to-sample fluctuations are very small. The analysis of the Ising model was performed 
to serve as a benchmark of the method. Recent high precision estimates, in the infinite volume limit, are as 
follows [fH^: T c = 4.511524(2), v = 0.63012(16). The results for K 3 = 0.8 and Ising confirm the 
effectiveness of the method, and suggest that the evaluation of irrelevant exponents is necessary in order for 
to obtain reliable results at higher precision. Notice the results from lattices (14, 10) : in this case a different 
transformation was used (the tie-breaker is not necessary). Our results for the Edwards-Anderson model 
are in agreement with previous estimates [|l3 14 1, obtained with finite-size scaling techniques on much 
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(24,8) 


4.4779 (2) 


0.6512(3) 




(24, 12) 


4.4767 (3) 


0.6439 (5) 




(24, 16) 


4.4769 (7) 


0.6395 (12) 


oo 


(12,8) 


4.5213 (2) 


0.6588 (9) 




(16,8) 


4.5178(1) 


0.6537 (7) 




(16,12) 


4.5149 (2) 


0.648 (2) 




(24,8) 


4.5149(1) 


0.6495 (5) 




(24, 12) 


4.5134(1) 


0.6441 (8) 




(24, 16) 


4.5128(1) 


0.6417 (14) 



Table 2: Estimates of the critical temperature (T c ) and the critical index v. K3 = oo denotes the Ising 
model (with periodic b. c). 

larger lattices and exploiting dedicated machines or supercomputers. 

Finally, we present in Table || the values of observables A\ , . . . . Aq at the estimated critical point. These 
data roughly characterize the fixed point on the lattice L' = 2. Remnant effects of irrelevant couplings 
affect the determination of T c , and therefore our extrapolation. Nevertheless, data points clearly point out 
the existence of different universality classes below and above the deconfinement transition. Within current 
precision, we can very well assume that model K3 = 0, 0.3 are in the same spin glass universality class, 
while K3 = 0.8 appears to be in the same class as Ising model. 

7 Conclusions 

A recently introduced jl8[ | (see however |||) Renormalization Group transformation suitable for disordered 
spin models is explained in detail and applied numerically in the three dimensional case. Universality 
classes depending on a frustration parameter have been characterized, pointing out the influence of de- 
confinement phase transition of gauge fields on the behaviour of the associated spin model, namely in the 
transition from glassy to ferromagnetic behavior. The study of universality classes relies on the direct detec- 
tion of fixed points, therefore on fundamental grounds, without relying on phenomenological comparison 
of critical indexes. Moreover, the RG analysis provided good estimates of the critical temperature, T c , and 
thermal index v, of the four models considered: K3 = 0, 0.3, 0.8 and Ising. As we have shown, the funda- 
mental determination of universality classes allows to consider vq and ^0.3 as independent estimates of the 
thermal index of the glassy fixed point. These results are in agreement with previous numerical estimates 
on the Edwards- Anderson model, obtained exploiting Finite-Size Scaling Ansatz Jl3| , |l4] , [16| |, For the small 
lattices and the reasonable computer time employed, this RG technique is to be considered very effective in 
order to compute critical properties of disordered systems. 



10 



K 3 




A 1 


A 2 


A A 


A 4 


^5 


A 6 





(12, 8) 
(16, 8) 
(16, 12) 


0.647 (6) 
0.650 (7) 
0.656 (18) 


0.580 (8) 
0.583 (9) 
0.589 (23) 


0.520 (7) 
0.523 (8) 
0.528 (20) 


0.653 (7) 
0.656 (8) 
0.660 (19) 


0.570 (8) 
0.573 (9) 
0.578 (21) 


1.143 (10) 
1.139 (11) 
1.133 (27) 


ft 1 


(16, 8) 
(16, 12) 


n fids ( 1 41 
0.670 (10) 
0.71 (3) 


ft Sid (171 

0.603 (12) 
0.65 (4) 


ft 51? { \ft\ 

0.543 (12) 
0.59 (4) 


ft (S4Q C 1 41 

0.673 (10) 
0.71 (3) 


ft SK7 Cl "SI 

0.592(11) 
0.64 (4) 


1 145 f?f)1 
1.114 (14) 
1.07 (4) 


0.8 


(12, 8) 
(16, 8) 
(16, 12) 
(24, 16) 


0.3429 (14) 
0.3505 (12) 
0.369 (5) 
0.367 (6) 


0.2920 (15) 
0.2990 (13) 
0.316(5) 
0.310(6) 


0.1840(11) 
0.1924(11) 
0.213 (4) 
0.215 (5) 


0.3969 (13) 
0.4032(11) 
0.419 (4) 
0.415 (5) 


0.2805 (13) 
0.2877 (11) 
0.305 (5) 
0.303 (5) 


1.5103 (21) 
1.5012(19) 
1.479 (7) 
1.489 (9) 


00 


(12, 8) 
(16, 8) 
(16, 12) 
(24, 16) 


0.3390 (4) 
0.3460 (3) 
0.3629 (6) 
0.3757 (8) 


0.2881 (4) 
0.2945 (4) 
0.3099 (8) 
0.3207 (9) 


0.1830 (3) 
0.1899 (3) 
0.2068 (6) 
0.2202 (7) 


0.3936 (4) 
0.3993 (3) 
0.4133 (6) 
0.4234 (8) 


0.2780 (4) 
0.2842 (3) 
0.2994 (6) 
0.3108 (8) 


1.5166 (6) 
1.5078 (5) 
1.4869 (11) 
1.4726(12) 



Table 3: Characterization of the fixed points: observables A\, . . . , A§ computed on the lattice Li at the 
estimated critical temperature T C (L\, L<i). 
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